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ABSTRACT 

Among 'empirical' methods of estimating oxygen abundances in extragalactic H II 
regions, the use of the ratio of nebular lines of [O in] and [N n], first introduced by 
Alloin et al. in 1979, is reappraised with modern calibration data and shown to have 
certain advantages over R 23 = ([O n] + [O III]) /H/3 and N2 = [N n] A6583/Ha, 
particularly when applied to star-forming galaxies at high redshifts. 
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1 INTRODUCTION 

Extragalactic H II regions are central to the study of chem- 
ical abundances in the universe (e.g. Garnett 2004). Accu- 
rate abundance measurements require the determination of 
the electron temperature, which in turn is usually obtained 
from the ratios of auroral to nebular line intensities, such as 
[O in] A4363/A5007. This is often referred to as the 'direct' 
T e method. A well-known difficulty, however, arises from the 
fact that, as the metallicity increases, the electron tempera- 
ture decreases (since the cooling is via metal lines) and the 
auroral lines eventually become too faint to measure. De- 
tailed modelling can overcome this problem. For example, 
Kewley & Dopita (2002) have recently used a combination 
of stellar population synthesis and photoionization models 
to develop a set of ionization parameter and abundance di- 
agnostics based only on the strong optical emission lines. 
Their 'optimal' method uses ratios of [N n], [O II], [O III], 
[S il], [S in] and the Balmer lines, that is the full complement 
of strong nebular lines accessible from the ground. They also 
recommended procedures for the derivation of abundances 
in cases where only a subset of these lines is available. 

The advent of 8-10 m class telescopes has made it pos- 
sible to extend observations of these lines to galaxies at high 
redshifts (e.g. Teplitz et al. 2000; Pettini et al. 2001; Kob- 
ulnicky et al. 2003; Lilly, Carollo, & Stockton 2003; Steidel 
et al. 2004). At high redshift, however, the study of nebular 
emission encounters new obstacles. The first and obvious dif- 
ficulty is that, given the much larger distances involved, the 
line fluxes are reduced dramatically. The problem is com- 
pounded by the fact that the lines are redshifted into the 
near-infrared wavelength region, where the sky background 
is orders of magnitude higher than in the optical. Second, 
at a given redshift only a subset of the strong lines will fall 
within an atmospheric window and thus be accessible from 
the ground. Even at particularly favourable redshifts, such 



as z ~ 2.3 which shifts [O n], [O ill] and H/3, and Ha and 
[N il] respectively into the middle range of the J, H and 
K bands, it is not possible to record all these lines in a 
single exposure. Rather, different spectrograph settings are 
normally required, and this can easily introduce additional 
errors in the relative flux calibrations. For all of these rea- 
sons, there are practical limitations to the accuracy with 
which emission line ratios can be measured in high redshift 
objects. 

There is thus an incentive to explore simpler abundance 
indicators based on only a few emission lines — preferably 
close in wavelength — which, while admittedly less accurate 
than the full treatments mentioned above, may still be ad- 
equate to characterise the chemical enrichment of distant 
galaxies. The idea of such indicators goes back to a fun- 
damental paper by Searle (1971) who, noting some mea- 
surements of nebular [O ill] lines in M33 by Aller (1942), 
himself found a systematic increase in [O ill] /H/3 and a cor- 
responding decrease in [N ll]/Ha with galactocentric dis- 
tance as a generic property of giant H II regions in late-type 
spiral galaxies. Searle pointed out that these trends were 
most likely due to radial abundance gradients, thereby rais- 
ing the possibility of using these line ratios more generally as 
indicators of abundance. Such indicators can be calibrated 
against more secure abundance determinations, based on au- 
roral line measurements or on detailed photoionization mod- 
els, to establish their typical accuracy; for this reason they 
are often referred to as 'empirical' methods. The first ex- 
ploration of such a method was by Jensen, Strom & Strom 
(1976) who considered the ratio [O ill] /H/3. The subsequent 
extension of this idea by Pagel et al. (1979), who included 
[O il], led to the most widely used abundance indicator, the 
R23 index which relates the abundance of oxygen to the ra- 
tio of ([O n] + [O ill]) to H/3. In the last 25 years there have 
been many efforts to refine the calibration of 7?23- The most 
successful are the calibration by McGaugh (1991) which is 
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Figure 1. Oxygen abundance against the N2 index in extragalactic H II regions. The symbols have the following meanings: Filled 
squares - (O/H) from photoionisation models by Diaz et al. (1991 — M 51) and Castellanos, Diaz, & Terlevich (2002a, b - NGC 925 and 
NGC 1637). Crosses - M 101 (Kennicutt et al. 2003). Open triangles - spiral and irregular galaxies (Castellanos et al. 2002a; Garnett 
et al. 1997; Gonzalez-Delgado et al. 1994, 1995; Kobulnicky & Skillman 1996; Mathis, Chu, & Peterson 1985; Pagel, Edmunds, & Smith 
1980; Pastoriza et al. 1993; Shaver et al. 1983 [30 Dor and NGC 346]; Skillman, Cote, & Miller 2003; Vflchez et al. 1988). Open squares 
- blue compact galaxies (Campbell, Terlevich, & Melnick 1986; Garnett 1990; Izotov, Thuan, & Lipovetsky 1994; Kniazev et al. 2000; 
Kunth & Joubert 1985; Pagel et al. 1992; Skillman & Kennicutt 1993; Skillman et al. 1994; Thuan, Izotov, U Lipovetsky 1995; van 
Zee 2000; Walsh & Roy 1989). Error bars in log (O/H) are only plotted when they are larger than the symbol to which they refer. The 
long-dash line is the best fitting linear relationship: 12 + log (O/H) = 8.90 + 0.57 X N2. The short-dash lines encompass 95% of the 
measurements and correspond to a range in log (O/H) = ±0.41 relative to this linear fit. The continuous line is a cubic function of the 
form 12 + log (O/H) = 9.37 + 2.03 X N2 + 1.26 X N2 2 + 0.32 X 7V2 3 which, however, gives only a slightly better fit to the data (95% 
of the data points are within ±0.38 of this line). The dot-dash horizontal line shows the solar oxygen abundance 12 + log (O/H) = 8.66 
(Allende-Prieto, Lambert, & Asplund 2001; Asplund et al. 2004). 



based on photoionisation models and the empirical one by 
Pilyugin (2003); both improve the accuracy by making use 
of the [O m]/[0 n] ratio as an ionisation parameter. Even so, 
systematic differences of up to 0.5 dex remain between the 
oxygen abundances derived from the R23 index and those 
deduced from the 'direct' T c method (e.g. Skillman, Cote, 
& Miller 2003; Kennicutt, Bresolin, & Garnett 2003). When 
compounded with the well-known double-valued behaviour 
of (O/H) vs. R23, this can lead to order of magnitude uncer- 
tainties in (O/H) in galaxies at z ~ 3 (Pettini et al. 2001). 

Partly to circumvent this latter problem, analogous 
S23 = ([S n] + [S m])/H/3 (Vflchez & Esteban 1996; Chris- 
tensen, Petersen, & Gammelgaard 1997; Diaz & Perez- 
Montero 2000) and & 3 4 = ([S 11] + [S m]± [S rv])/H/3 
(Oey et al. 2002) indices have been proposed. Although these 
sulphur-line based indices turn over at higher metallicities 
than J?23, the sulphur lines are weaker than their oxygen 
counterparts and their red and infrared rest wavelengths 
place them beyond the reach of ground-based observations 
at redshifts 2 £ 1.5. 

More promising are methods which involve [N 11] and 
Ha. In this paper we reconsider the N2 index (N2 = 
log{[N n] A6583/Hq}) recently discussed by Denicolo, Ter- 
levich, & Terlevich (2002), and show how its accuracy can 



be further improved in the high abundance regime by com- 
parison with [O in] A5007/H/3. 



2 THE N2 INDEX 

Following up on earlier work by Storchi-Bergmann, Calzetti, 
& Kinney (1994) and by Raimann et al. (2000), Denicolo et 
al. (2002) have focussed attention again on the N2 index, 
pointing out its usefulness in the search for low metallic- 
ity galaxies. The [N n]/Ha ratio is highly sensitive to the 
'metallicity', as measured by the oxygen abundance (O/H), 
through a combination of two effects. As (O/H) decreases 
below solar, there is a tendency for the ionization to in- 
crease (either from hardness of the ionizing spectrum or 
from the ionization parameter, or both), decreasing the ra- 
tio [N n]/[N ill], and, furthermore, the (N/O) ratio itself de- 
creases at the high-abundance end because of the secondary 
nature of nitrogen (e.g. Henry, Edmunds, & Koppen 2000; 
Pettini et al. 2002a). From the ground, Hcv and [N 11] can 
be followed all the way to redshift z ~ 2.5 (the limit of 
the K band). With thousands of galaxies now known at 
z > 1 (Madgwick et al. 2003; Steidel et al. 2004; Abraham 
et al., in preparation), the N2 index offers the means to de- 
termine metallicities in a wholesale manner over look-back 
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Figure 2. Oxygen abundance against the 03N2 index in extragalactic H II regions. The symbols have the same meaning as in Figure 1. 
The long-dash line is the best fitting linear relationship, 12 + log (O/H) = 8.73 - 0.32 X 03^2, valid when 03^2 < 1.9 The short-dash 
lines encompass 95% of the measurements which satisfy this condition and correspond to a range in log (O/H) = ±0.25. The dot-dash 
horizontal line shows the solar oxygen abundance 12 + log (O/H) = 8.66 (Allende-Prieto, Lambert & Asplund 2001; Asplund et al. 2004). 



times which span most of the age of the universe. It is thus 
worthwhile reconsidering its accuracy. 

Denicolo et al. (2002) compiled an extensive sample of 
nearby extragalactic H II regions which could be used to cal- 
ibrate A^2 vs. (O/H). We have revised their database by: (i) 
including only H II regions with values of (O/H) determined 
either via the T e method or with detailed photoionization 
modelling, and excluding the untabulated estimates for H II 
regions from the catalogue of Terlevich et al. (1991), because 
for many of these (O/H) had been estimated from 'empir- 
ical' (i.e. R23 or 523) indices; (ii) updating it with recent 
measurements by Kennicutt et al. (2003) for H n regions 
in M 101 and by Skillman et al. (2003) for dwarf irregular 
galaxies in the Sculptor Group. 

The total sample, shown in Figure 1, consists of 137 ex- 
tragalactic H 11 regions with well determined values of (O/H) 
and N2. In all but six cases (shown with filled symbols) the 
oxygen abundance was determined with the T e method. Ref- 
erences to the original works are given in the Figure caption. 
It can be readily appreciated from Figure 1 that, while A?" 2 
does indeed increase monotonically with log (O/H), a linear 
relationship is only an approximate representation of the 
data, which rather tend to lie along an S-shaped curve. In 
particular, when (O/H) ~ 8.2±0.2, the value of N2 can span 
one order of magnitude, from A^2 ~ —1.8 to N2 ~ —0.8, 
whereas at solar metallicity and beyond [N 11] tends to sat- 
urate (Baldwin, Phillips & Terlevich 1981) because it comes 
to dominate the cooling (Kewley & Dopita 2002). Forcing a 
linear fit to the sample, 1 we obtain a line of best fit according 
to the relation: 



12 + log (O/H) = 8.90 + 0.57 x N2 (1) 

shown with the long-dash line in Figure 1. Both the slope 
and the intercept are lower than the values of the fit by 
Denicolo et al. (2002) who proposed 12 + log (O/H) = 9.12 + 

0. 73 x N2. While the formal statistical errors on the slope 
and intercept are small (0.03 and 0.04 respectively), of more 
interest is the dispersion of the points about the line of best 
fit in Figure 1. Specifically, 95% (68%) of the measurements 
of log (O/H) lie within ±0.41 (±0.18) of the line defined by 
eq. (1); the 2a limits are shown as short-dash lines in Figure 

1. Only a marginally better fit to the data is provided by a 
third order polynomial of the form: 

12±log (O/H) = 9.37+2.03xAf2+1.26xAf2 2 +0.32xAf2 3 (2) 

(valid in the range -2.5 < N2 < -0.3), with 95% (68%) of 
the measurements being within ±0.38 (±0.18) of the values 
given by eq. (2). This cubic fit is indicated by the continuous 
line in Figure 1. We conclude that with the N2 calibrator it 
is possible to estimate the abundance of oxygen to within a 
factor of ~ 2.5 at the 95% confidence level. This accuracy is 
comparable to that of the R23 method. 



3 THE OW2 INDEX 

We now consider to what degree the accuracy in the determi- 
nation of (O/H) can be improved by considering two ratios, 
[N II]/Hq and [O m]/H,fl. Alloin et al. (1979) were the first 
to introduce the quantity OSN2 = log{([0 ill] A5007/H/3)/ 



1 We performed an unweighted least squares fit to the data be- 
cause in compilations such as this the error estimates are highly 



heterogeneous. However, a weighted least squares fit gives very 
similar results. 



4 M. Pettini and B.E.J. Pagel 



([N n] A6583/Hq)}, 2 but since then the 03N2 index has 
been comparatively neglected in nebular abundance stud- 
ies. There is now a sufficient body of high quality data to 
reassess its merits. We expect the inclusion of [O ill] to be 
most useful in the high metallicity regime where [N n] sat- 
urates but the strength of [O ill] continues to decrease with 
increasing metallicity. 

Figure 2 shows how 03iV2 varies with (O/H) for the 137 
extragalactic H n regions in our sample. Clearly this method 
is of little use when 03JV2 ^ 2, but at lower values there 
appears to be a relatively tight, linear and steep relationship 
between OSN2 and log (O/H). A least squares linear fit to 
the data in the range —1 < 03N2 < 1.9 yields the relation: 

12 + log (O/H) = 8.73 - 0.32 x 03iV2 (3) 

which is shown by the long-dash line in Figure 2, with 95% 
(68%) of the measurements within 0.25 (0.14) dex of the best 
fit line; the short-dash lines in Figure 2 show the 2a limits. 
The caveats are that the statistics are somewhat limited, 
with only 65 out of the 137 H n regions in our sample sat- 
isfying the condition 03JV2 < 1.9, and that the relationship 
defined by eq. (3) relies heavily on the four data points with 
12 + log (O/H) ~ 9.0 deduced from photoionisation mod- 
els rather than the 'direct' T c method. Nevertheless, in the 
high metallicity regime the OSN2 method appears promis- 
ing and it would be very worthwhile extending the database 
of nearby H II region measurements to improve its statistics. 



4 CONCLUSIONS 

In this note we have reassessed the usefulness of 'empirical' 
methods based on the [N ll]/Ha ratio for determining the 
oxygen abundance in H II regions, particularly with an eye 
to their application to the analysis of star-forming galaxies 
at high redshift. We have revised the extensive compilation 
prepared by Denicolo et al. (2002) to include only H n re- 
gions where the oxygen abundance is believed to be known 
reliably, either because the electron temperature has been 
determined with high precision, or through detailed mod- 
elling of the ionisation of the nebula. Our main conclusions 
are as follows. 

1. The relationship between N2 = log ([N II] A6583/Ha) 
and log (O/H) is only approximately linear and yields (with 
the revised data set considered here) estimates of (O/H) 
which are only accurate to within ~ 0.4 dex at the 95% 
confidence level. This level of precision is comparable to that 
of the more widely used R23 method. 

2. At moderate to high metallicities, greater than ~ 1/4 
solar, the oxygen abundance can be deduced to within a 
factor of ~ 0.25 dex (again at the 95% confidence level) if 
03N2 = log{([0 in] A5007/H/3)/([N 11] A6583/Ha)} is lower 
than ~ 1.9. The OSN2 index is particularly useful at solar 
and super-solar metallicities, where N2 saturates. 

When applied to high redshift galaxies, the N2 and 
03iV2 indices offer several advantages compared to the 
more familiar R23 method. First, the N-based indices have 

2 This definition is slightly different from the original one pro- 
posed by Alloin et al. (1979) who included both [O III] doublet 
lines in the numerator of the first ratio. 



a monotonic behaviour with (O/H), overcoming the ambi- 
guities introduced by the double- valued character of i?23- 
Second, both N2 and 03N2 rely on ratios of emission lines 
which are close in wavelength, so that corrections for red- 
dening, and accurate flux calibration of the spectra, are not 
necessary. Third, both indices are highly sensitive to the 
oxygen abundance and are thus particularly well suited to 
the analysis of data of only moderate signal-to-noise ratio. 
Specifically, the slope of 0.57 of the relationship between 
log (O/H) and N2 implies that one only needs to measure 
the [N ii]/Hq ratio with an accuracy of about 25%. Bet- 
ter precision would not improve the estimate of (O/H) since 
the random error would then be of minor importance rela- 
tive to the systematic uncertainty of ±0.4 dex in the method. 
Similar considerations apply to the 03JV2 index which is in- 
herently more accurate in the high metallicity regime and 
varies even more steeply with log (O/H). 

The high sensitivity of OSN2 to the oxygen abundance 
also helps to mitigate another complication which may arise 
preferentially at high redshift. The galaxies accessible to 
detailed spectroscopic studies at redshifts z — 1 — 3 sup- 
port very high rates of star formation, one to two orders of 
magnitude higher than those seen in the local universe (e.g. 
Steidel et al. 2003, 2004). There are also mounting indica- 
tions that, even at these early epochs, such galaxies have 
already reached near-solar metallicities (e.g. Pettini 2004). 
When galactic chemical evolution proceeds at such a fast 
pace, the production of primary nitrogen from intermedi- 
ate and low mass stars can lag behind that of oxygen from 
Type II supernovae, as found by Pettini et al. (2002b) in 
the z = 2.7276 Lyman break galaxy MS1512-cB58 (see also 
Matteucci & Pipino 2002). How far (N/O) is 'out of equi- 
librium' depends on the metallicity, the star formation rate, 
and the yet poorly constrained timescale for the release of 
primary nitrogen (e.g. Meynet & Pettini 2003). However, 
even in cases where N is over-deficient relative to O by fac- 
tors of 2-3, as in MS1512-cB58, the slope of -0.32 in the 
relation log (O/H) vs. OSN2 limits the error in log (O/H) 
to ^ 0.15 (given that ([N 11] /Ha) varies linearly or slower 
with the nitrogen abundance in this regime). 

Finally, it is encouraging that abundance calibrators 
developed from measurements in individual extragalactic 
H 11 regions seem to be fairly robust when applied to the 
integrated light of galaxies (Kobulnicky, Kennicutt, & 
Pizagno 1999; Moustakas & Kennicutt, in preparation). 
Thus, we propose that, with the aid of the 03JV2 index, it 
should be possible to measure the oxygen abundances of 
star forming galaxies at z > 1 to better than a factor of two 
with only a relatively modest investment in observing time. 

We are grateful to Chris Akerman, Mike Irwin and 
Samantha Rix for help with the production of the figures, 
and to Mike Dopita, Lisa Kewley and Chuck Steidel for valu- 
able comments. 
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